Direct observation of room temperature high-energy resonant excitonic effects in 

graphene 
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Using a combination of ultraviolet- vacuum ultraviolet reflectivity and spectroscopic ellipsometry, 
we observe a resonant exciton at an unusually high energy of 6.3eV in epitaxial graphene. 
Surprisingly, the resonant exciton occurs at room temperature and for a very large number of 
graphene layers N^75, thus suggesting a poor screening in graphene. The optical conductivity 
((Ti) of resonant exciton scales linearly with number of graphene layer (up to at least 8 layers) 
implying quantum character of electrons in graphene. Furthermore, a prominent excitation at 
5.4eV, which is a mixture of interband transitions from tt to tt* at the M point and a n plasmonic 
excitation, is observed. In contrast, for graphite the resonant exciton is not observable but strong 
interband transitions are seen instead. Supported by theoretical calculations, for A'^ < 28 the cri is 
dominated by the resonant exciton, while for A'' > 28 it is a mixture between exitonic and interband 
transitions. The latter is characteristic for graphite, indicating a crossover in the electronic 
structure. Our study shows that important elementary excitations in graphene occur at high 
binding energies and elucidate the differences in the way electrons interact in graphene and graphite. 



PACS numbers: 

Graphene, a one layer thick carbon honeycomb struc- 
ture, has recently attracted a lot of attention due 
its exotic quantum properties[T]. In addition to the 
well-known electronic properties, such as ballistic elec- 
tron transport [TH3j: quantum Hall effect [1], tunable 
band gapO E], and physics driven by many-body 
interactions [7], graphene also displays very interesting 
optical properties. For instance, due to its low-energy 
excitation graphene becomes highly transparent in the 
visible spectral range and its infrared optical conductiv- 
ity can be tuned using a gate voltage [3]. These particular 
properties of graphene potentially result in exciting op- 
toelectronic applications. 

Recent theoretical studies based on the ab-initio GW 
and Bethe-Salpeter equation (BSE) approach by Yang et 
al.|10j and Trevisanutto et al.fTT] have predicted the ex- 
istence and have highlighted the importance of resonant 
excitonic effects in the optical absorption of graphene. 
However, there is disagreement as to the origin and po- 
sition of the exciton. In Ref. [TD] the calculations were 
done up to 7 eV and the resonant exciton predicted to oc- 
cur at 4.6 eV due to electron-hole interaction in tt*- and 
TT-band at the M point. While in Ref. [TT] the calcula- 
tions were done at a much higher energy of 22 eV and the 
resonant exciton predicted to appear at 8.3 eV due the 



background single particle continuum of dipole forbidden 
transition at the the F point. Despite their disagreement, 
both have agreed that resonant exciton plays important 
role for elementary excitations in graphene and thus its 
understanding is crucial. 

A direct way to probe resonant exciton in graphene 
is to measure its complex dielectric response in wide en- 
ergy range pT]. Despite reports on the optical proper- 
ties of graphene over the visible energy ranged |H1 HH- 
there has been no reliable optical conductivity data 
at higher energies. Therefore, it is crucial to study the 
optical properties of graphene in an unprecedented high 
energy range not only to provide a direct evidence of the 
resonant excitonic effects per se, but also for our under- 
standing of electronic structure of graphene. 

Herewith, we report the optical conductivity (cti) of 
epitaxial graphene on 6H-SiC(0001)/buffer layer sub- 
strate in an unprecedented wide photon energy range 
from 0.5 to 35 eV using a combination of spectro- 
scopic ellipsometry and UV-VUV reflectrometry |17| . 
It has been recently shown that the combination of 
these two techniques enables a stabilized Kramers-Kronig 
transformation ^1 7j which is crucial to resolve precisely 
the dielectric function. We study the evolution of cti as 
a function of the number of layers TV (henceforth CTi^at) 



2 




5 10 15 20 25 30 35 
Photon energy (eV) 



10 15 20 25 30 35 
Photon energy (eV) 



FIG. 1: (color online) Room temperature experimental re- 
sults of (a) reflectivity (b) loss function, Im (e~^). The inset 
of (a) shows the experimental geometry while the inset of (b) 
shows (e^^) on an expanded scale from 4.5 to 6.5 eV and from 
5 to 10 eV. 



including graphite. 

Epitaxial graphene samples on 6H-SiC(0001) were pre- 
pared by Si desorption process in Si-flux as described in 
previous studies [T51 [T^]. This process was found to lead 
to multilayer graphene with ordered stacking TJ*. The 
thickness was determined using scanning tunneling mi- 
croscopy and Raman spectroscopy [191 [20]. For reference, 
we have used a cleaved highly ordered pyrolytic graphite. 
The reproducibility was checked by fabricating and re- 
peating individual measurements on three different sam- 
ples that were grown under identical conditions at least 3 
different locations on each sample yielding reproducible 
results. The details of the samples are described in the 
supplementary section. 

Figure [ija) shows room temperature reflectivity data 
for A'^-dcpendent epitaxial graphene, graphite as well as 
the substrate. Interestingly, one can see that the reflec- 
tivity of epitaxial graphene is very much dependent on A'^ 
and very much distinguished from that of graphite. The 
reflectivity has rich and distinct structures especially in 
the range from 5 to 8 eV for A'^ = 1, 2, and 8 while 
additional structures occur around 14 to 20 eV for 
= 22, and 75. On the other hand, the reflectivity of 
graphite has shown strong structures below 0.5 eV, ~ 5 
eV, and around 14 to 20 eV which are similar to published 
result Further rigorous discussion will be achieved 
by analysing the cti. (Details of analysis of reflectivity, 
complex dielectric function and optical conductivity are 
shown in the Supplementary.) 

The ai^N shows striking results (Figs. [2] (a) and (b)). 
For ~ 1, the <Ji i is dominated by two well-defined 
peaks, very pronounced at 6.3 eV (peak B) and less pro- 
nounced at 5.4 eV (peak A), with almost equal full-width 
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FIG. 2: (color online) The optical conductivity (ui) of (a) 
substrate and graphene (A*' = 1,2,8) and (b) graphene (A'' = 
22,75) and graphite show three peaks at 5.4 eV (label A), 
6.3 eV (label B) and 14.1 eV (label C), respectively. Inset 1 
shows the comparison of (Ti,i between experimental data and 
theoretical calculations. Insets 2 and 3 show the cri on an 
expanded scale at various energy ranges, (c) The value of cti 
at A, B and C as a function of N. (d) The partial spectral 
weight (W) as a function of N. 



at half maximum (FWHM) of ~0.8 eV. These structures, 
especially the peak B, are considerably sharp for this high 
energy range. 

Now, we shall focus our discussion on peak B as 
it shows an interesting dependence on A^ (Fig- [2jc)). 
The ai^N{B) increases rapidly towards A^=8, while its 
FWHM and peak position remain independent of A^. For 
A^ > 8, (7i_Ar(B) decreases while the peak position shifts 
towards higher energy (6.7 eV and 11 eV for N = 22 and 
N = 75, respectively). For graphite, peak B disappears 
and our result is similar to the published data[2lj. 

To find out the origin of peak B, we have di- 
rectly compared our experimental data with theoretical 
calculations|10 [ fTT j . As shown in inset 1 of Fig.[2|^a), the 
line shape and the CTi i(B) between experimental data 
and theoretical calculations from Ref. [TT] are surpris- 
ingly very similar. This comparison shows a decisive ev- 
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idence of the high energy resonant exciton in epitaxial 
graphene. The resonant exciton arises from dipole tran- 
sitions of the single-particle continuum. However, the 
observed resonant exciton peak occurs red shifted com- 
pared to the calculations. Our calculations show that 
the red shifted exciton could result from film-substrate 
interactions, which are not included in the previous cal- 
culations (see below for the details). Thus, the origin of 
peak B is the high energy resonant exciton as predicted. 

Furthermore, our detail study rules out the interband 
transitions and plasmonic excitations as origin of peak 
B. Based on Density Functional Theory (DFT), we have 
calculated optical conductivity for iV = 1 and 2, and 
found that the peak B does not originate from interband 
transitions (see Fig. 3a and discussion below) . Secondly, 
one may argue that we have to consider contributions 
from a plasmonic excitation. To address this issue, we 
have studied in detail an energy loss function Im(e~^) 
which can reveal collective excitations such as plasmonic 
excitations [221 [53]. As shown in Fig. 1(b), there is no 
feature that can be attributed to a plasmonic excita- 
tion. For = 1, Im(e~^) is dominated by a strong 
structure at ^10 eV and a weak structure at 5.5 eV. 
These structures are from tt + cr and tt in-plane plasmon 
modes, respectively |23i . Based on the symmetry, the 5.5 
eV plasmon peak structure is only visible for light po- 
larization parallel to c-axis (E||c). In our measurement, 
the E was mixed between E||c and E_Lc with main con- 
tribution from E_Lc, thus the observed tt plasmon is very 
weak. On the other hand, tt + cr plasmon mode is very 
strong in our geometry which is consistent with EELS 
measurements [251 [2^ . As N increases, the tt + a plas- 
mon mode shifts toward higher energy and gets broader, 
while the structure at 5.5 eV is nearly A-dependent. This 
blueshift may be due to strong effects of the interlayer 
Coulomb interaction on the total T)lasmon[22|. 

The exciton energy of 6.3 eV suggests an untypically 
high binding energy which can also be made responsi- 
ble for the presence of the exciton up to room tempera- 
ture. If one considers the angular-resolved photoemission 
data, the distance of a band at F point to Fermi level is 
around 4.5 eV[3S]. By assuming the distance of a* band 
to Fermi level to be similar to that of a band one will 
end up with 9 eV as the separation between cr and cr* 
bands. Based on these values the binding energy of the 
resonant exciton is 2.7 eV which is 108 times larger than 
the thermal energy. This binding energy is significantly 
larger as compared to that of ZnO which is around 60 
meV (2.4 times thermal energy) [27] and is very clearly 
seen at room temperature while the binding energy of 
the excition for GaAs is about 3.4 meV and can only be 
seen at low temperatures [58J. Such a very high binding 
energy exciton observed on graphene is indeed remark- 
able and a direct consequence of its low dimensionality 
in conjunction with its intriguing electronic properties. 

We next discuss peak A and its origin. As shown in 
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FIG. 3: (color online) (a) comparison between experimen- 
tal data and calculated ai. Upper inset : cti spectra of 
graphite. Lower inset : cti spectra of the substrate, (b) 
The background single-particle continuum transitions (A) at 
the r point as a function of nominal charge transfer (rie) from 
substrate to graphene for various graphene layer (A). Red ar- 
row shows the observed charge transfer in our sample, (c) A 
proposed model of optical absorption of exciton in graphene 
and graphene on substrate as function of A'^. RE stands for 
the resonant exciton and (Aexc) is excitation energy. 



<^i,n{A) increases monotonically as a function 
of A^. Interestingly, for A^ > 8 as well as for graphite 
the position of peak A shifts towards lower energies as 
opposed to peak B, while the line shape is rather sym- 
metric and the FWHM is broader for larger A^. Based on 
electronic band structure calculations f IT, peak A is the 
result of interband transitions from the tt to tt* bands at 
M-point where the van Hove singularity occurs. Further- 
more, from analysis of energy loss function (Fig. [ijb)) 
the Im(£^^) shows a weak 7r-plasmon contribution. Thus 
one can conclude that the peak A is a mixture between 
interband transitions and a plasmonic excitation. 

Another notable observation is a broad structure at 
14.1 eV (peak C). The peak C is absent for A^ = 1 and 
CTi_Ar(C) increases dramatically for N > 2 (Fig. [2jc)). 
Interestingly, for graphite the cti is dominated by peak 
C with symmetric lineshape. Based on electronic band 
structure calculations [HI [HI |25] this structure is inter- 
band transitions from cr to a* bands at F point. Peak C 
clearly has an intimate relationship with optical conduc- 
tivity in graphite while peak B is a unique characteristic 
of graphene. 

We next discussed partial spectral weight integral (W) 
because it describes the effective number of electrons 
excited by photons of respective energy. The cri is re- 
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stricted by the f-sum rule: ai{E)dE — where 
n is the electron density, e is elementary charge and m* 
is effective electron mass. Hence, one can extract the 
W = J^^ ai{E)dE, for various energy ranges. Due to the 
f-sum rule and charge conservation, the W is constant, 
and thus one can study the spectral weight transfer and 
reveal interactions as well as the effect of N in the broad 
energy range of 0.5 to 35 eV. Figure ^d) shows W for 
the different energy regions: 0.5 to 3.2 eV (region I, Wj), 
3.2 to 5.8 eV (region II, Wn), 5.8 to 8.9 eV (region III, 
Win), and 8.9 to 35.0 eV (region IV, Wiv)- 

The Wi is mainly governed by the transition around 
Dirac cone, i.e., tt to tt* around K point in Brillouin 
zone which is consistent with previous publications [71 [T5] . 
While the W77, W///, and Wiv show the main contribu- 
tion from the peak A, B, and C, respectively. Interest- 
ingly, Wi and Wji show almost similar trend in which 
the W increases monotonically as function of N, while 
the Win and Wiv show completely different fashion. 
For < 8 the Wm increases as N increases while for 
> 8 Win decreases as A^ increases. On the other 
hand, Wiv increases rapidly for N > 8. In graphite, 
Wiv is maximum while Wm obtains its smallest value. 
We find a crossover between Wm and Wiv at A^ ^ 28 
(see Fig. [gd)). Thus, we propose that for N < 28 
the optical conductivity of epitaxial multilayer graphene 
is dominated by high energy resonant excitonic effects, 
while for A^ > 28, the optical conductivity is dominated 
by interband transitions. 

To gain further insight, we have calculated the opti- 
cal conductivity (ai^caic) of graphene, graphite, and sub- 
strate using Density Functional Theory (DFT) . We have 
studied the role of interband transitions and have com- 
pared them with experimental results (see Fig. |3][a)). 
Noting that our calculations do not include electron-hole 
(e — h) interactions, thus {(Ji^caic) is mainly driven by 
interband transitions and serves as reference for the non- 
correlated case. (The detail of calculation is given in 
supplementary section.) Based on the DFT calculation, 
io'i,caic) mimics reasonable well the (iTi) for graphite as 
well as for 6II-SiC(0001)/buffer layer. For graphite, the 
peak at 4 eV is dominated by transitions from tt n* 
bands while the peak at 14.1 eV is transitions from 
cr — > (7* bands. This is consistent with previous theo- 
retical study [25 . 

On the other hand, the DFT calculations for graphene 
show that its optical conductivity cannot be explained 
with interband transitions. As shown in Figure |3][a), 
the calculated optical conductivity for graphene shows a 
completely different result than the experimental data 
and thus our attempt to mimic cti for A^ = 1 and 2 
failed. Calculations based on DFT show that the (Ti, caic 
for graphene (A^ = 1, 2) above 3 eV is very similar to 
graphite. This is in contrast to our experimental results. 
This further supports that the (Ti^n is mainly driven by 
strong e — h interactions which lead to high energy reso- 



nant exciton and form unusual electronic band structure 
while the optical conductivity of graphite is mainly driven 
by interband transitions. It is shown that for A^ = 1, once 
one turns on the e — h interactions, the peak around 4 
eV and peak around 14.1 eV (inset of Fig. |3]^a)) vanish 
resulting new and very strong peak in between [TT]. 

To find out the origin of the redshift of the excitonic ex- 
citation, we have calculated the variation of the a bands 
and the background single-particle continuum transitions 
(A) as function of charge transfer (rie) and strain. We 
find out that while the cr bands are nearly independent 
from these two effects, the background single-particle 
continuum transitions are very much dependending on 
the charge transfer rif, from substrate to the graphene 
layer (Fig. [sjjb)). Based on our DFT calculations, the 
charge transfer from substrate to graphene rig is about 
^0.07e (per graphene unit) and thus the background 
single-particle continuum transitions reduce by ~1.4 eV 
compared to that of for n,, = 0. The strain in graphene 
due to the lattice mismatch with the substrate reduce 
the background single-particle continuum transitions by 
0.6 eV. These two effects altogether decrease the single- 
particle continuum transitions to about 2.0 eV in total. 
Thus, one may expect to see the exciton at lower energy 
of ~6.3 eV. Our result may suggest that energy excita- 
tion of exciton depends on the rig and the strain. 

Interestingly, our calculation shows that A does not 
dependent on A^ significantly (Fig. [3][b)). This is in fact 
consistent with our experimental result. In contrast, the 
theoretical predictions in Ref.[TTj show an energy shift 
between the exciton in graphene and bilayer graphene 
which is attributed to enhanced screening. Our exper- 
imental results show, however, that the position of the 
exciton does not depend on the A^ value till about to 
N = 22. This can be reconciled with the Ref.[TT] only 
if one assumes that screening effects are significantly re- 
duced. Thus, by comparing the experimental results and 
theoretical calculations, one may conclude that screen- 
ing effects in graphene are much weaker than one would 
expect. Altogether our results suggest that graphene be- 
longs to the class of strongly correlated electron systems. 

Finally, based on our experimental results of high en- 
ergy optical conductivity, DFT calculations, and recent 
ab initio GW-BSE calculations [TT], we have proposed 
the following phase diagram (Fig. [2]^c-d)) and a model 
for optical absorption which is proportional to cri (Fig. 
|3]jc)). For free standing graphene of N=l, the optical 
absorption is dominated by high energy resonant exciton 
effects which occur around 8.3 eV. However, in our case 
due to the charge transfer and lattice mismatch between 
graphene and the substrate, the resonant exciton reduces 
to 6.3 eV. For A^ < 28, the optical absorption shows very 
strong resonant excitonic-like structure and weak inter- 
band transitions. While for 28 < A^ < 75 the optical 
absorption shows a mixture between exciton and inter- 
band transitions (Figs. [2]^d) and |3|c)). For A^ > 75 as 
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well as for graphite, the optical absorption is dominated 
by interband transitions. 

In summary, we have observed a high energy room tem- 
perature stable resonant exciton at ~6.3eV with a large 
binding energy of ~2.7 eV in the optical conductivity 
of multilayer epitaxial graphene revealing strong collec- 
tive e — h interactions. The resonant exciton persists for 
very large N and thus dominates the electronic properties 
owing to the poor screening in graphene. Furthermore, 
the mixture of interband transition from n to it* at the 
M point and weak tt plasmonic excitation gives rise to 
peak in the optical conductivity at 5.4 eV. These find- 
ings demonstrate the importance of high energy optical 
conductivity and have strong implications on the under- 
standing of the electronic structure of epitaxial graphene. 
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